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Abstract 

We show the possibihty of achieving experimentally a Takens-Bogdanov bifurcation 
for the nonlinear interaction of two transverse modes (/ = ±1) in a CO2 laser. 
The system has a basic 0(2) symmetry which is perturbed by some symmetry- 
breaking effects that still preserve the Z2 symmetry. The pattern dynamics near 
this codimension two bifurcation under such symmetries is described. This dynamics 
changes drastically when the laser properties are modified. 
PACS number(s):42.65.Pc, 02.20.-a, 42.55.Lt 
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1 Introduction 



Pattern formation in physical systems is an area of active research. It has been 
recently pointed that lasers can provide a "test bench" for these studies, as it is 
possible to have many active transverse modes which, through nonlinear interaction, 
give a complicated spatiotemporal dynamics [|l|. Recent evidence shows that some 
interesting phenomena such as defects, vortices, chaotic alternancy, etc., can be 
found in laser systems . Unfortunately, in most cases these works do not enlighten 
the mechanisms behind the transition from simple solutions to the spatiotemporal 
uncorrelated ones. A first step in understanding this complex dynamics is to study 
the nonlinear interaction among a few modes which bifurcate from the zero solution 
|]. This article strives in this direction. 

It has been reported that some qualitative features of the transverse patterns of 
intensity observed in CO2 lasers can be explained in the frame of the theory of bifur- 
cations in the presence of symmetries [0, ^. The laser tube imposes the 0(2) symme- 
try (rotations and reflections). The mode amplitude equations having this symmetry 
predict that the stable modes with non-zero angular momentum bifurcating from 
zero should be travelling waves, in contradiction with experiments. Anisotropics in 
the laser parameters (pumping, losses or disalignments in the setup, etc.) break the 
0(2) symmetry. An agreement between theory and experiment can be reached by 
introducing a symmetry-breaking term preserving the Z2 symmetry. This model 
succesfuUy predicts the stability of the standing waves widely observed in the laser 
transverse section. The prediction of secondary solutions that arise from the stand- 
ing waves is another achievement of the model with "imperfect" symmetry. These 
solutions are either a mixture of a standing wave and travelling waves, or modulated 
waves that appear from a Hopf bifurcation of the standing waves 

The conditions for a OO2 laser displaying those solutions that bifurcate simul- 
taneously from the standing wave (codimension two point) are analyzed here. The 
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dynamics near such bifurcation and the intensity patterns resuhing from this sit- 
uation are investigated. Work in this direction has been done in traveUing-wave 
convection in bynary mixtures. In this case, the presence of distant sidewalls in 
systems that are translation- invariant break the 0(2) symmetry and the simplest 
possible symmetry-breaking effects are discussed when the system undergoes a Hopf 
bifurcation 0. 

This paper is organized as follows. In section 2 we discuss the properties of the 
physical system and we present the equations describing the dynamics of the active 
primary modes in this system. A linear study of these modes is performed in order to 
identify the possible secondary solutions that might bifurcate from them. In section 3 
we analyse the conditions under which a nonlinear interaction between the secondary 
solutions is possible, and a normal form reduction of the equations is carried out. 
After computing the coefficients of this normal form in terms of parameter values 
for a CO2 laser, we determine the two possible scenarios . We discuss the kind 
of patterns that could be observed under these conditions in section 4. Section 5 
sets out some conclusions and provides some guidelines for further experimental 
observation. 



2 The Model and Primary Bifurcations 



The physical system is a CO2 laser in a Fabry-Perot cavity. The active medium is 
contained in a cylindrical tube, with a perfectly reflecting plane mirror at one end 
and a curved mirror with partial reflectivity at the other end. Physically the effec- 
tive curvature of this mirror can be modified by inserting a passive optical device. 
Moreover, we are interested in the interaction among modes with nonzero angular 
momentum. This can only be achieved experimentally by placing an intracavity iris 
which inhibits the Gaussian mode 

We therefore consider that the electric field can be expressed in terms of the 
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modes arising from a Hopf bifurcation from the zero solution: 



E = P,{r)L,{z)e'^\z^e''' + z^e-'^'), (1) 

where Pi (r) and Li (z) account for the radial and longitudinal dependence of the 
bifurcating modes, respectively, and 6 is the angular variable. The integer / denotes 
the angular momentum and uj is the temporal frequency of the modes. The complex 
amplitudes zi and Z2 are governed by the equations 



Xzi - A{zizl + 2z2zl)zi + ez2 



4 = \z2 - Ai2zizl + Z2zl)z2 + ezi 



(2) 
(3) 



which are obtained by substituting the electric field E into the Maxwell-Bloch 
equations and truncating to third order |^. The complex coefficients A and A 
can be expressed in terms of the atomic inversion and decay rates, the detuning 
and a convolution between the pumping profile and the spatial part of the modes. 
(The explicit form of these coefficients can be in the Appendix). Note that the 
complex parameter e (which carries information about the asymmetries in the laser 
parameters, such as anisotropics in the pumping or in the Brewster windows [Q) 
accounts for the breaking of the 0(2) symmetry (see Appendix). As the equations 
remain invariant under the operation zi <-> they preserve a Z2 symmetry. 

We substitute Zi = pje*'^% e = p^e'^'^' and \ = p + iu in equations (2-3). After 
scaling by A^ = Re{A) we arrive to the following system: 



P'l = I^Pi - (Pi + 2p2)Pi + PeP2 cos(5 + (p,) (4) 

P'2 = PP2 - (2pi + P2)P2 + PePl C0S(5 - 0^) (5) 

S' = a(pi -P2) -Pe(Pi/P2sin(5-0,) +p2/pisin(5 + 0e)) (6) 
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and a fourth equation for the evolution of the phase 0i which is uncoupled from these 
three equations. The variable 5 = 02 — 0i is the phase difference between the two 
modes. Notice that its dynamics is nontrivial because the broken symmetry term 
couples this phase difference with the mode amplitudes. (For simplicity we use the 
same notation for the rescaled parameters). The information about the curvature, 
detuning, etc. contained in remains in the coefficient a = A^/A". 

The primary bifurcations from zero in equations (4-6) are four stationary so- 
lutions: two standing waves SWo{pi = p2,S = 0), SWt,{pi = P2,S = tt) and 
two solutions which are a mixture between traveling waves and standing waves 
TWi{pi > p2,6 > 0), TW2{pi < p2,S < 0). These solutions give four stationary 
patterns their stability depending on the parameter settings. By an appropriate 
selection of the parameters in this situation a good qualitative agreement with a 
previous experiment has been obtained [Q]. 



3 Takens-Bogdanov Bifurcation 



In this paper the attention is focussed in investigating the secondary bifurcations 
from the standing wave and on the new dynamics arising from the nonlinear in- 
teraction between those secondary solutions. Before proceeding with the stability 
analysis it is convenient to perform the following change of variables: 



pi = 5cos(0/2) (7) 
P2 = 5sin(0/2). 

because the B direction is decoupled from the (0, 6) plane in this representation 
(Fig.l). The standing wave solutions are simply SW^y. 



2 u±p^cos0. 
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A linear stability analysis around these solutions leads to the following eigenval- 
ues. In the B direction the eigenvalue is simply — 2(// ± pecos(0e)), where the ± 
correspond to the SWo^^^ solutions respectively. The eigenvalues corresponding to 
the {(j), S) plane are obtained from the following matrix: 

|;U =1= cos0e ±2p^sm(f)^ 

Pe COS 4>e) T 2pe sin 0^ T'^Pe COS 0^ 

Now we determine the conditions for the appearance of a codimension-two (CT) 
point. 

3.1 Codimension-Two Point 

The codimension of a bifurcation is the smallest dimension of a parameter space wich 
contains the bifurcation in a persistent way. In our case the CT point is determined 
by a degenerated double zero eigenvalue in the matrix (8). This results from the 
interaction between a Pitchfork bifurcation and a Hopf bifurcation. (The Pitchfork 
bifurcation gives rise to two new solutions {TW{^2)- The solution resulting from the 
Hopf bifurcation is the so-called modulated wave (MW)). This CT point is obtained 
when the parameters verify: 

p — ±5peCOS0g (9) 

a-i = -tan(2(/),) (10) 

For a given a, and as p > (the two modes were born), these relations (9- 
10) lead to four solutions as sketched in Fig. 2. Solutions I and IV are associated 
with the SWo (sign + in eq. 9): (I = 0°; 71/ = 0° + 37r/2), while II and III (77 = 
0°+7r/2; 777 = (/>°+7r) correspond to SWj, (sign - in eq. 9). As cqs.(4-6) are invariant 
under the transformation (0^; pi{t), p2{t), S{t)) — >• (0g -|- tt; pi{t), p2{t), d{t) + tt) {Z2- 
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residual symmetry) the local dynamics around I and III is equivalent. The same 
argument can be applied to solutions II and IV. 



For the sake of simplicity we analyze only the local behavior around instability 
I(and IV) in the next subsection. (Calculations for cases II (and III) follow the same 
procedure). 

3.2 Normal Form Reduction 

Let us recall some general features of a dynamical system linearized around a sta- 
tionary solution. When some eigenvalues have zero real part, the flow near this fixed 
point can be quite complicated. The linear space spanded by the states correspond- 
ing to this null real part eigenvalues is known as center eigenspace. The invariant 
manifold tangent to the center eigenspace constitutes the center manifold. The local 
dynamics 'transverse' to this manifold is relatively simple, since it is controlled by 
the 'fast' variables of the flow. The asymptotic behavior of the flow develops on the 
center manifold. The family of these behaviors that arise in the vicinity of this fixed 
point when the parameters are slightly varied is called the unfolding of the bifurca- 
tion. The 'simplest' set of equations that reproduces generically all these behaviors 
is called the normal form of the bifurcation. 

The unfolding of the bifurcation for solution I (dynamics near SWq) is described 
by the normal form: 



that was first studied by Takens and Bogdanov [§, 3. it can show in fact two 
different behaviors depending the choice of the sign =f. 

Now we show that eqs.(4-6) reduce to eqs. (11-12) around SWq if the conditions 
(9-10) for solution I hold. 



x' = y 



(11) 



y' = ax + by ^ — x^y 



(12) 
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Let us note that in this case 0° G (0,7r/2) for any given a (eqs.9-10). We 

2 2 

perform a last change of variables: {B,(f),S) —>■ {B,u,S) with u — — — cos0. 

The equations can be rewritten: 

B' = iiB -]-{3-u'^)B^ + p,B{l-u'^)^/^ cos (peCos 5 (13) 

u' = m(1 - + 2p£(l - ?i^)^/^(sin(;6csin(5 - ?icos</),cos(5) (14) 

5' — aB^u . (cos 0e sin (5 + K sin cos (5) (15) 

yl — ii 

To obtain the unfolding of the bifurcation we consider for small variations of the 
parameters: 

— — 5 cos (f)^ + q 

Pe 

The fast variable B is 'enslaved' by the other two variables and can be adiabatically 
eliminated after assuming B' = 0. Then the central manifold is two-dimensional and 
can be expressed as a function of (u,5). As usual in a local stability analysis, we will 
keepup to third order terms. To obtain the normal form (eq. (11-12)) the following 
change of variables (u,5) ^ (x = Kto^ - S),y ^ ^^^^ + '^)) introduced. 
This allows to obtain the right linear part of the normal form. We can rearrange 
the equations by means of another near-identity change on the new variables {x, y). 
After these changes the normal form of the bifurcation reads as: 



where 



x' ^ y (16) 

y' = nix + H2y + Cx^ + Dx^y (17) 



4 l + ata(7(0O) atag^cj)^,) + 8atag{(j)^^) + 7 
- 3^ cos(0O) ^ (18) 



eJ 



(19) 
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D 



C 



a 




(22) 



(20) 



(21) 



Notice that the coefficients in this normal form mainly depend on the value of a. 
As D is always positive the unfolding of the bifurcation depends mainly on C wich 
is a function mainly of a. The unfolding is qualitatively different when C changes 
from negative to positive, i.e., for the conditions: 



(The unfoldings corresponding to C > and C < will we described in detail 
in the next section). In the appendix the dependence of a on the laser properties is 
calculated for different transverse modes. This calculations show that the transition 
point (eq. (23)) can be reached in a CO2 laser (Fig. 3). 

4 Pattern Dynamics in the Neighborhood of the 
Bifurcation 

In section III we obtained the conditions for the two posible kinds of T-B bifurcation 
that can take place in the transverse section of a CO2 laser. Depending on the sign 
of parameter C two qualitatively different scenarios are possible. Now we describe 
the two different dynamics in some detail: 

a) ae (-00, 0.258) = (f), e (0, 0.912): 

This corresponds to the case C > (D < 0). The unfolding near the bifurcation 
is presented in Figure 4a. From this unfolding the following general features are read. 
First, the solutios SWq (Fig. 5) as well as the two TW (created after a Pichtfork 




(23) 
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bifurcation from the SWo) are unstable for nearly all the parameter space. The 
system almost ever evolves to some of the travelling solutions, TWi = {pi,p2,5) 
or TW2 = {p2, Pi,—5). So, one of the conjugate patterns in figure 5 would be 
observed when we tune the parameter values around the bifurcation point. Second, 
the analysis of the full equations (4-6) shows that far away from the bifurcation point 
a limit cycle appears as a global bifurcation from a heteroclinic connection between 
these patterns {TWi^2)- This solution will look like a periodic alternancy between 
the two TWi^2 due to the critical slowing down that takes place in the neighborhood 
of these points (Fig. 5). 

h) q; e (0.258, 00) = 0e e (0.912, f): 

In this case C < (L> < 0) and there is a strong change in the pattern evolution 
(Fig. 4b). Near the bifurcation point the SWq, weak oscillations TWI 2 ^ SWo 
(quasi-stationary patterns), or a periodic alternancy TW{ ^ SWo ^ TW2 will be 
observed. The latter corresponds to a limit cycle that grows in a global bifurcation 
from a homoclinic connection of the SWo- When we go far away in the parameter 
space a periodic alternancy between the two TWi^2 is found again. 

Notice that the frequency associated to the oscillation between the patterns 
(Hopf bifurcation from 5*1^0) is of order ujI^^ ~| A'-i |~ Pe | 5 |- So it is much slower 
than the temporal scale associated to the modes. 

5 Conclusions 

In this article, a model for the evolution of transverse modes (/ = ±1) in a CO2 
laser has been studied. Due to the unavoidable anisotropies in the laser setup, some 
terms that break the natural 0{2) symmetry must be included. This model is quite 
succesful in predicting the stability of the primary solutions SW found in recent 
experiments. The possible secondary solutions have been determined. 
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The symmetry breaking term allows the possibility of a co dimension- two point 
bifurcating from the SW. This point appears when secondary solutions from a 
Pitchfork and a Hopf bifurcation interact simultaneously at the same point (Takens- 
Bogdanov bifurcation). A normal form reduction procedure allowed to capture the 
main dynamical properties of the system near that bifurcation. The calculations of 
the normal form coefficients for different realistic values in a CO2 laser showed that 
two possible dynamical scenarios are possible. 

The intensity patterns that correspond to the dynamical interaction between 
these modes (/ = ±1) have been described for these two scenarios. Near the bifurca- 
tion one can get: stationary patterns (standing waves SW or travelling waves TW) 
or oscillations between the SW and the new conjugate patterns TW'. Far away 
from the point bifurcation a global bifurcation of a heteroclinic connection leads to 
a periodic alternancy between TW patterns. These patterns look quasi-stationary 
due to a critical slowing down near the TW solutions. 

Remarkably this " rich" dynamics cannot be observed without a symmetry- breaking 
term. It is important to stress that the interpretation of complicated pattern dy- 
namics beyond the particular problem studied here could lie on simple symmetry 
considerations. 
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7 Appendix 



It is clear that the analysis presented in this work strongly relies on the terms 
multipliying e in Eqs. (2-3). An asymmetry in the pumping profile gives rise to 
this linear coupling between the equations for the mode amplitudes zi and Z2. The 
dependence of e in terms of laser parameters was derived in reference . The result 
was: 

m 



where is the parameter governing the asymmetry in the pumping profile {Kpumping - 
K{r) + 2Kf, cos(2^)), (3 is the rate of decay of the atomic polarization and fli is the 
slow temporal frequency of the empty cavity mode. 

In the following we compute the value of a = ^ for the primary bifurcating 
modes with angular momentum / = ±1. It is shown that the transition point 
a = 0.258 between the two possible T-B bifurcations can be obtained in a CO2 
laser. 

The spatial coordinates of the problem are (r, 6, z) where (r, 9) corresponds to 
the transversal section of the cavity and z to the longitudinal direction {z G [—L, 0]). 
The new coordinates {^,6,z) are introduced 0: 



,2R„ 



ksis + z 



2\-l 



1)1/2 



(25) 
(26) 



where s is the effective curvature, L(~ Im) the cavity length, the curvature 
radius of the spherical mirror and k = is the wave number. 

The field is expanded (following [0) in terms of empty cavity modes: 



(27) 
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This expansion is introduced in the Maxwell-Bloch equations (a brief review is done 
in IQ), and projected on each mode. The temporal evolution for the cavity mode 
amplitudes Za is then: 

Z'o, = KZa + Maf^^pZ^Z^Zp + h.o.t. (28) 

where the expressions of the coefficients A^, Ma^^p are given in 
The mode = (a+, a~) has the functional form: 



a, = (R^iOe-^^'e^''', R,iOe"^-e-''-') (29) 
The coefficient (prnu is given by: 



(f)^ = Q^z-p^tan{-) (30) 

s 

with 

~ n^vr + arctan(-) — (5 (31) 



where {p^, l^, n^) are three integers that characterize the mode: = 2nr + + 1 
is essentially the total tranverse energy of the laser beam [rir is the radial quantum 
number), is the angular momentum around the z axis and is associated to 
the longitudinal behaviour of the mode. -R/i(0 is the radial dependence and 6 = 
k mod 2tt is the detuning of the cavity. 

Making the operations indicated in 0, we obtain the general form: 

T 2 

M^^^p = F{n^,n,,np)— I J J Ki^R^iOR^iOMOMO x (33) 

[4 cos(</)„ + (pu - 4>,i - 4>i3) + 2 cos(0a - <j)u - 4>f, + (f)p)]^d^dzd6 
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where the pumping profile K[^) = K is considered constant along the whole beam 
width and x gives account for the cavity losses. 

After applying these procedures to our two mode interaction model: 



z[ = Xzi - A{zizl + 2z2zl)zi + ezi (34) 

4 = )^Z2 — A{2zizl + Z2zl)z2 + ez2 (35) 

we can identify A — — Mim that, according to expression (31) gives: 

^ = _M„„ = cte.^^/6fl.md« (36) 

These calculations have been performed for the primary bifurcating modes — 
(p,l,n): 



6i^(2,l,0) ^ nm-ee'^' (37) 
ci^(4,l,0) ^ Rl(^)^(2-efee-^' (38) 
61^(2,1,1) ^ nm-ee-^' (39) 



In general the coefficient a — A^A^ depends on the cavity curvature s and on 
the detuning 5. We plot the dependency of a on s and 5 for the primary modes in 
figure 3. a decreases when the curvature (or detuning) increases and the transition 
value a = 0.258 can be reached for a sufficiently high curvature s. 
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Figure Captions 



1. Sketch of the change of variables (7) performed in the equations (4-6). In 
these variables the dynamics of interest takes place near the {6, 4>) plane (center 
manifold) . 

2. The four possible solutions of the equation system (9-10) (under the condition 
> 0) where a Takens-Bogdanov bifurcation takes place. They are shown in the 

angular (p^ space. The dynamics around all them is equivalent. 

3. (a) Plot of the value of a as the curvature s of the cavity is varied for three 
different modes: bi = (2,1,0), ci = (4,1,0) and b[ = (2,1,1). The detuning 5 is 
fixed: 5 = 0. 

(b) Plot of a on s for the mode bi — (2, 1, 0) when the detuning takes the values: 
5 = 1, 5 = 2 and 5 = 3. 

(c) Plot of a on 5 for three modes with different longitudinal behavior: bi — (2, 1, 0), 
b[ = (2, 1, 1) and b'l = (2, 1, 2). (The curvature is fixed: s = 20). 

The transition explained in relation (23) is obtained for a — 0.258. (See appendix 
for the notation). 

4. Unfolding diagrams for the two kind of Takens-Bogdanov bifurcation (see 
the normal form 16-17) that can be achieved: squares= TW[ 2 (or TWi^2), circles= 
SWq. The outer figures of the diagram are obtained simulating the total number of 
equations (4-6) far from the instability. 

(a) For a < 0.258 where C > 0, L> < 0, 

(b) For a > 0.258 where C < 0, D < 0. 

5. Patterns that can be observed: TWi (or TW[ if it is close to SWq), SWq and 
TW2 (or TW^2)- The arrows represent the periodic alternancy TW[ <-> SWq <-> TW2 
that takes place for C < 0. It is also posible another periodic alternancy TWl ^ 
TW^ or TWi ^ TW2 for C> 0. 
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